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Abstract 

A new method for extracting neutron densities from intermediate energy elastic proton-nucleus 
scattering observables uses a global Dirac phenomenological (DP) approach based on the Relativis- 
tic Impulse Approximation (RIA). Data sets for ^'^Ca, ^^Ca and ^''^Pb in the energy range from 
500 MeV to 1040 MeV are considered. The global fits are successful in reproducing the data and in 
predicting data sets not included in the analysis. Using this global approach, energy independent 
neutron densities are obtained. The vector point proton density distribution, is determined 
from the empirical charge density after unfolding the proton form factor. The other densities, /?^, 
P^, Ps, are parametrized. 

This work provides energy independent values for the RMS neutron radius, Rn and the neutron 
skin thickness, in contrast to the energy dependent values obtained by previous studies. In 
addition, the results presented in paper show that the expected rms neutron radius and skin 
thickness for ^''Ca is accurately reproduced. The values of Rn and 5„ obtained from the global fits 
that we consider to be the most reliable are given as follows: for '^'^Ca, 3.314 > i?„ > 3.310 fm and 
-0.063 > Sn> -0.067 fm; for ^^Ca, 3.459 > Rn > 3.413 fm and 0.102 > Sn > 0.056 fm; and for 
208p|^ 5.550 > Rn > 5.522 fm and 0.111 > Sn > 0.083 fm. These values are in reasonable agreement 
with nonrelativistic Skyrme Hartree-Fock models and with relativistic Hartree-Bogoliubov models 
with density-dependent meson-nucleon couplings. The results from the global fits for '^^Ca and 
^^^Pb are generally not in agreement with the usual relativistic mean-field models. 

PACS numbers: 25.80.Dj, 24.10.Jv, 24.10.Ht, 21.60.-n 
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I. INTRODUCTION 



Determination of the proton and neutron densities, their root-mean-square radii, Rp and 
Rn, and the neutron skin thickness, S'„ = i?„ - Rp, are critical to understanding many of 
the bulk properties of matter [l|, Q, y, Q ■ Horowitz et al. have pointed out that there are 
substantial disagreement between theoretical values of 5^ jsf. Furnstahl's recent analysis 
of neutron radii in the framework of mean-field models shows that relativistic mean-field 
models overestimate the values of S'n Q. Parity violation electron scattering may provide 
the experimental data to resolve the differences in the theoretical values 0], but is also useful 
to have alternative methods of obtaining neutron densities. Reliable neutron densities are 
needed for atomic parity violation experiments jHHSIS Hi; the analysis of antiprotonic 
atoms in understanding the surface crust of neutron stars [l^, and in extrapolation 
to proton- rich or neutron-rich nuclei that is important in nuclear astrophysics 0]. In this 
work, we revisit the analysis of medium energy proton-nucleus elastic scattering data with 
the goal of obtaining reliable, energy independent neutron densities and the values of R^ 
and Sn- The analysis of elastic electron scattering which has resulted in reliable ground 
state charge densities has been a guiding light for our work [l^ . 



For a number of years we have used the relativistic im pul se approximation (RIA) in the 
analysis of proton- nucleus elastic and inelastic scattering |l6l IitI isl lol . 2^ 21 , 2^, l^] and 
the RlA-KDP(Kemmer-Duffin-Petiau) H Hi, HI for meson-nucleus elastic scattering [27|, 



|28, . „2S. .30. .,31l] . These approaches produce relativistic optical potentials which result in good 
agreement with medium-energy scattering observables. The input to these calculations are 
the relativistic densities from Quantum Hadrodynamics (QHD) js^, S^] and the elementary 
NN amplitudes from Arndt et al. [s^]. In recent work we use the modern EFT densities 

in H mill 113 • 

The seminal analysis of proton-nucleus elastic scattering data done by Ray and Hoff- 
mann used both the RIA and the non-relativistic KMT approach in their fits to get the 
observables from 300 MeV to 1040 MeV |1^. Unfortunately neither approach produced 
energy-independent neutron densities. In addition, some of the values of Sn for ^^Ca and 
208pb were negative, in contradiction to all nuclear structure calculations. Shlomo and Schaf- 
fer used the results from an analyses of 1 GeV proton elastic scattering from ^''Ca and ^^Ca 
to obtain the skin thickness for ^°Ca, ^^Ca, '^^Ca and ^®Ca, see Table 2 in J39]. Starodubsky 
and Hintz extracted the neutron densities from elastic proton scattering of ^''^'^^'''^^'^Pb at 650 
MeV and obtained Sn of (0.20 ± 0.04) fm for ^ospi^ However, the energy-independence 
of the neutron densities in the work of Ref. |1^ or Ref. was not addressed. Recently 
Karataglidis et al. have calculated proton and neutron elastic scattering from ^°^Pb and 
■^^Ca targets at three energies 40 MeV, 65 MeV and 200 MeV. They used a model based on 
coordinate space nonlocal optical potentials using a full folding NN interactions with various 
Skyrme model ground-state densities i^]. For ^osp]^ they found that the SKM* model gave 



the best agreement with proton and neutron elastic scattering data at 40 MeV, 65 MeV and 
200 MeV. Based on this the authors suggest that Sn for ^°Ca is ~ -0.05 fm and ^osp^ jg ^ 
0.17 fm 41]. 

In this paper, a new analysis of proton-nucleus elastic scattering is used to obtain the 
neutron density. This work is motivated by our considerable experience in obtaining high 
quality global proton-nucleus optical potentials from 20 MeV to 1040 MeV |3]. The new 
method meshes the global approach with the RIA and proves to be successful in obtaining 
energy independent neutron densities. The starting point is the RIA in its simplest form. 



2 



which for spin-zero nuclei includes only scalar, vector and tensor terms. The tensor term 
is very small and is excluded; as was done in the RIA analysis done by Ray and Hoffmann 

We find that substantial progress in extracting the neutron densities from proton-nucleus 
elastic scattering is made by using a global approach focusing on the energy region where 
the RIA is capable of reproducing experiment very well. We have obtained values of Sn 
for ^''Ca, ^^Ca, and ^°^Pb which agree with nonrelativistic Skyrme models 0, 0, li^ and 
relativistic Hartree-Bogoliubov model extended to include density dependent meson-nucleon 
couplings Our results for '^^Ca and ^osp]^ g^^.^ generally not in agreement with relativistic 
mean-field model, see Ref. [6] and references therein. 

The organization of this paper is as follows. Section |n] describes the global method 
to obtain the neutron density. Section IIIII discusses the results and the sensitivity of the 
extracted neutron density, i?„ and Sn to the input used in the fitting procedure. The values 
of Sn and the neutron RMS radii i?„ for "^^Ca, '^^Ca, and ^osp]^ f^j. ^-j^g various tests of the 
input to the model are given in this section. The summary and conclusion is in Sect. IIVI 



II. RIA GLOBAL METHOD FOR EXTRACTING THE NEUTRON DENSITY 

The RIA nuclear reaction formalism is used as the basis of global fits to medium-energy 
proton nucleus elastic scattering data. The input to the RIA consists of the Arndt NN 
amplitudes and the point proton density, which is fixed from the charge distribution 
obtained from electron nucleus scattering. The neutron vector density, and the scalar proton 
and neutron densities are parametrized, resulting in good fits of p + A elastic scattering 
data between 500 MeV and 1040 MeV. Using the RIA as the basis for the global fits is a new 
approach and our results shows that it is a valid method for extracting neutron densities, 
Rn and Sn- 

In the global approach used in this work the form of the RIA vector and scalar and tensor 
optical potentials are given by 

Usir) = E f Wdq^Mqr)F^{q)pi{q), (1) 

where the fourier transforms of the density form factors are 

pi(g) = I dVe^^^-V.(r'), (4) 
pliq) = [ dVe^^^->i(r'), (5) 
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m = I rfVe^^^-^lp^ (r'). (6) 

The subscripts s,v,t refer to Lorentz scalar, vector (time-like), and tensor quantities. The 
superscript n and p refer to neutrons and protons, F{q) are the invariant NN amplitudes, 
and R{q) is the kinematical factor required to obtain the invariant NN amplitude in the 
Breit frame |38i]. The value of the upper limit on the momentum transfer, q, is determined 
by the available on-shell NN data. Many studies conducted have demonstrated that higher 
order corrections to this first order RIA-Dirac optical model approach are negligible at the 
energies being studied in this paper. See Ray et al. for more details of the RIA-Dirac 
calculations used for elastic scattering observables [l^ . 

Elastic p + A data between 500 MeV and 1040 MeV for ^^Ca, ^^Ca and ^o^Pb form the 
data set. The quality of the fits are good and the predictions of data not in the data set is 
used to verify the procedure. However, unlike the usual RIA, where we have generally used 
scalar, vector and tensor densities from the results of relativistic effective field theory (EFT) 
calculations, we use the RIA as a basis for extracting the vector neutron density as well as 
the two scalar densities using a global fitting procedure. 

Next we describe the treatment of the vector and scalar densities. The point proton 
density is fixed by using the results from electron scattering as follows. 



where, 



P'.ir) = ^ / d'qe-'^-p'M (7) 



_ J d'qe^^^-%{r) 
P'M) - — • (8) 



Here pc{r) is the experimental charge density, G{q) is the proton form factor. The point 
proton density, p^(r) is normalized to Z and the neutron density p"(r) is normalized to N. 

The next step is to consider the model that will be used. After several years of testing 
a number of different forms used to parametrize the neutron vector and scalar densities, 
we have chosen to model the vector neutron density and the scalar proton and neutron 
densities using the cosh form (COSH) used in many of our global fits j4^. We have found 
the cosh form parametrization produced more stable results for R — n and Sn than the three 
parameter fermi (3PF) or the sum of gaussian (SOG) parametrizations. This stability of 
the cosh was understood from the results from many global fits using the all three forms. 
In global fits we used four or five different momentum transfer ranges. For example, if five 
momentum transfer ranges are used the fits are done from 0.0 fm~^ to 1.5 fm~^ to 0.0 fm~^ 
to 3.5 fm~^ in steps of 0.5 fm~^. In future work we will consider other parametrizations. 

The form of the cosh model in this work the superscript V stands for volume and S for 
surface terms, 

|cosh[it/aJ + cosh[r/aJ — 2} 

5 {cosh[R/a] - l}{cosh[r/a] - 1} 

^ ^ ' ' ^ {cosh[i?/a]+cosh[r/a]-2}2 ' ^ 
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The vector neutron density is 



p^^{r, R,, a,) = p^(r, R,, a^) - p^(r). (11) 

The density p^{r, Rb, ab) is given by 

p^{r,Rb,ab) oc {r, Rb,ab) + af^ {r, Rb,ab), (12) 

and p^{r,Rb,ab) is normahzed to A. There are two geometric parameters, Rb and ab, and 
the parameter, a. 

The scalar proton and neutron densities are, 

p?(r, i??, <) oc fir, Rl, <) + (3f{r, <) (13) 

and 

p:(r, R:, <) oc /^(r, i?:, <) + 7/^ (r, i^^, <). (14) 

Each of these densities contains three parameters, two geometry parameters, i?^, i?^, a^, 
and the parameters /3 and 7. The 10th parameter, Pio, searched is given by J (i^rp^ = 
Pio Z. A similar 11th parameter given by j d^rp^ = Pu N could also have been searched. 

However, the searches were more stable if the ratio j^^^y^, i-e., the ratio of was fixed 
by ratio of the volume integrals per particle of the scalar proton and neutron densities from 
any one of the EFT models given by Rusnak et al. [36.. .;37,] . Three different EFT densities 
are used, two are point coupling models, VA3 and FZ4, and one is a meson model, MA4. 
The parameter Pio is not sensitive to the EFT model chosen and the final difference and the 
scalar proton and neutron densities are not sensitive to the EFT model chosen even though 
the values of Sn for N 7^ Z in these models differ widely. For example, values of Sn for ^°^Pb 
are 0.332 fm for VA3, 0.259 fm for MA4 and 0.160 fm for FZ4. 

As mentioned above we found that the model used in this paper was the best of many 
different models we tried. We will investigate other models in future work using this global 
approach. 



III. FITTING PROCEDURE AND RESULTS 



energies 
MeV EjJ 
11 




It has long been known that using a global approach has been very useful in obtaining NN 
amplitudes and the NN phase shifts. This is one of the reasons for using a global approach 
when the data set used is large and usually correlated. In this work we use elastic p + A 
data between 500 MeV and 1040 MeV for ^^Ca, ^^Ca and ^ospb. For ^^Ca there are five 
497.5 MeV 0, 613 MeV HJ, 650 MeV 0^97.5 MeV ^ and 1040 
For 208pb there are five energie^ 497.5 MeV BH^, 613 MeV jij, 650 MeV 
55L797.5 MeV |56U57|, and 1040 Me VjsaL However for "^^Ca only three gies; 497.5 
MeV j45|, 797.5 MeV j59|, and 1044 MeV [5l|,|52| are available. In order to make predictions 
of data not included in the data set we remove one energy. For '^^Ca and ^°^Pb the 650 MeV 
data has been excluded, and for ^'^Ca 1044 MeV data has been excluded. The quality of 
the global fits to the data are good and the predictions for data not in the global data sets 
verify the procedure. 

In previous Dirac phenomenological (DP) global fitting, the data sets were cut at 100 
degrees in the center of mass or at momentum transfer q at 3.0 fm~^, whichever came first. 
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In this work we do global fits using a variety of momentum transfer values as discussed 
below. In addition, as is discussed later, we also test the sensitivity of the input to the 
model. 

In the global fitting we have used five momentum transfer ranges from 0.0 fm~^ to 1.5 
fm^^ to 0.0 fm~^ to 3.5 fm~^ in steps of 0.5 fm~^. Or four momentum transfer ranges from 
0.0 fm~^ to 2.0 fm~^ to 0.0 fm~^ to 3.5 fm~^ again in steps of 0.5 fm~^. The values of Sn and 
Rn are expected to change with the value of the momentum transfer range used as the data 
sets are changed. We find that the values of Sn and Rn for a given momentum transfer range 
(we use five momentum transfer ranges or four momentum transfer ranges) are clustered in 
a reasonably well constrained set of values for Sn and we call these stable results. The 
momentum transfer range for 0.0 fm~^ to 1.5 fm~^ is almost always the outlier, which is 
understandable as the data set is small and does not have the diffraction structure needed 
when fitting proton scattering from nuclear targets. Of course in any such global approach 
the data sets are usually correlated. This is certainly true as the data sets sets with different 
momentum transfer ranges do overlap. Thus a statistical analysis to obtain a mean value 
and a standard deviation can not generally be used. However, we can determine the range 
of the values of 5*^ and Rn for each fit from every momentum transfer range, and obtain the 
values for Sn and i?„ which could be shown as a scatter plot, as is usually done when data 
sets used are correlated. We have obtained the range of values for Sn and Rn for every test 
considered, while these ranges can not be interpreted as resulting from a statistical analysis 
we found that doing such an analysis gives a useful guide to present the results. 

There are a number of features in the fitting procedures which could produce changes in 
the neutron density, Rn and Sn- As mentioned above, two ranges of momentum transfer 
cuts are used. In the following sensitivity tests we use both ranges. Generally the smaller 
range produces the smallest range of values for Rn and Sn so in most of the figures we show 
the larger range values for Rn and Sn which gives the most conservative results. 

We test the results obtained when the following features of the model used in the fitting 
procedures are changed. First the three different EFT models used in fixing the 10th param- 



eter are considered [36|,|37[. To investigate the effect of the G{q) used in obtaining the point 
proton density we use two different forms they are identified as Gi{q) from [60] and G2{q) 
from [g^I • The set chosen for the Arndt NN amplitudes is input to the fitting procedure, 
and we use sets FAOO and SM86 to find the sensitivity of the fits to this input [3^ . The data 
sets included in the fits have been changed to see if the results change significantly. Finally, 
we use the three different charge density models used in the analysis of electron-nucleus 
elastic scattering data; the three parameter fermi (3PF), the sum of gaussian (SOG) and 
the fourier-bessel (FB) jH^I- Thus each target we have done global fits which have used 
effective field theory densities MA4, FZ4 and VA3, and for each three charge distributions 
are used. In addition two different form factors Gi{q) and (^2(9) and two sets of Arndt 
amplitudes, set FAOO and set SM86. Figures 1-5 show the results of these tests of the global 
fitting procedure. 

The results for the Rn and Sn are given in three tables, found in the Appendix; Table I 
for 40Ca, Table II for ^^Ca and Table III for ^ospb. 

There are three general cases we use in testing the global fits: Case 1 uses Arndt am- 
plitudes set FAOO and the Gi{q) form factor; Case 2 uses Arndt amplitudes set FAOO and 
the G2{q) form factor; and Case 3 uses Arndt amplitudes set SM86 and the Gi{q) form 
factor. In every case, we obtain good global fits for MA4, FZ4 and VA3 for every charge 
distribution used; SOG, 3PF, and FB. As discussed above the momentum cuts ranges used 
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are 0.0 fm~^ to 1.5 fm~^ to 0.0 fm~^ to 3.5 fm~^ in steps of 0.5 fm~^ or we remove the 0.0 
fm~^ to 1.5 fm^^ set. For each case the one with five momentum transfer ranges or four 
momentum transfer ranges we calculate the mean and standard deviation for i?„ and Sn for 
every global fit, this gives us an average over all EFT models for a given charge distribution. 
This is denoted by AVE eft in the tables. This is not to be taken as a statistical error, it is 
however a convenient way to show the range of the values in a consistent way. We do this 
rather than remove every outlier. Then for each case we show the combined values of the 
range for i?„ and Sn for that case, this gives the average over all charge distributions as well 
as over all EFT models. These results are denoted as AVEcd and AVEeft in the tables. 
Finally, we calculate the ranges of -R„ and Sn for all three cases combined for each target. 

The last calculation, which combines all of the values of Rn and Sn for every test made, 
gives us the most conservative results. From these final test cases, shown at the bottom of 
the tables in the appendix, we get the following values using the five momentum transfer 
ranges for ^°Ca, 3.350 > i?„ > 3.300 fm and -0.008 > 5„ > -0.080 fm; for ^^Ca, 3.505 > 
Rn > 3.421 fm and 0.148 > 5„ > 0.058 fm; for ^ospb, 5.589 > Rn > 5.513 fm and 
0.156 > 5*^ > 0.076 fm. As mentioned earlier if we use the four momentum transfer ranges 
results are generally more clustered and we obtain: for ^°Ca, 3.350 > Rn > 3.304 fm and 
-0.006 > Sn > -0.080 fm; for ^^Ca, 3.485 > i?„ > 3.417 fm and 0.129 > 5„ > 0.053 fm; 
and for ^ospb, 5.561 > Rn > 5.513 fm and 0.130 > 5„ > 0.074 fm. The values for i?„ and 
Sn for both of the final case results are quite similar. 

However, the authors found that the global fits using the Arndt NN amplitude set FAOO 
and the form factor Gi{q) gave the most stable results i.e. the values did not vary very 
much, see Tables I, II and III in the Appendix. This is especially true for using the SOG 
charge distribution. In this case the values of Rn and Sn using the five momentum transfer 
ranges are: for ^^Ca, 3.318 > Rn > 3.308 fm and -0.059 > Sn > -0.069 fm; for ^^Ca, 
3.498 > Rn> 3.412 fm and 0.141 > Sn > 0.055 fm; and for ^ospb, 5.602 > i?„ > 5.512 fm 
and 0.164 > Sn > 0.074 fm. The values for i?„ and Sn using the four momentum transfer 
ranges are: for ^^Ca, 3.314 > i?„ > 3.310 fm and -0.063 > S'„ > -0.067 fm; for ^^Ca, 
3.459 > Rn> 3.413 fm and 0.102 > S'„ > 0.056 fm; and for ^ospb, 5.550 > Rn > 5.522 fm 
and 0.111 > ^„ > 0.083 fm. 

The results using the form factor G2{q) are quite similar as well as the case using form 
factor Gi{q) but with the NN set SM86. These values as well as the values usingthe obtained 
from the final case values are in agreement with nonrelativistic Skyrme models ^, 4^, 4^, the 
relativistic Hartree-Bogoliubov model extended to include density dependent meson-nucleon 
couplings 01 and the recent analysis of antiproton atoms jl2| . Our results are generally not 
in agreement with relativistic mean-field models. 

Next we discuss the results of each test of the model separately. Figure [T] shows the 
effect of the two different momentum transfer ranges as well as the difference in the three 
EFT models used to fix parameter 10. The filled boxes (the smaller momentum transfer 
range) and diamonds (the larger momentum transfer range) are the values from the three 
different EFT models used. The charge density used in this case was the SOG for all 
targets, the proton form factor Gi{q), and the set FAOO Arndt NN amplitudes. This figure 
shows that the results from the two different momentum ranges overlap and that the EFT 
model used does not affect the results. The figure also shows several theoretical results from 
Refs. BiBBHHHH |4l|. The comparison between theory and range of the values 



from the fits shows that all theoretical Sn values agree with the results of the global fits for 
'^^Ca. However, the theoretical Sn values for ^^Ca and ^°^Pb are somewhat larger than the 
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global fits. 

Next we consider the results when two different proton form factors, Gi{q) from Ref. |60j 
and G2{q) from Ref. joH, are used in obtaining the vector point proton density for a given 
charge distribution model. The five momentum transfer ranges, the values for the AVEeft 
models for the SOG model charge distribution and set FAOO Arndt NN amplitudes are used 
and shown in Figure |21 

Investigating the sensitivity to the set of Arndt NN amplitudes used is done by comparing 
sets FAOO and SM86. The larger momentum transfer range, the averaged EFT models for 
the SOG model charge distribution, and the proton form factor Gi{q) are used and shown 
in Figure El 

In order to check the sensitivity due to the data sets included in the fit we have done 
fits using only two data sets (497.5 MeV and 797.5 MeV) for ^°Ca and ^ospb. The results 
agree very well with the same case using four data sets as is shown in Fig. EJ The values 
of Sn for the two data set case are for ^°Ca, —0.054 > Sn > —0.066 fm; and for ^osp^^ 
0.166 > Sn > 0.088 fm. These results agree very well with the values for the four data set 
case ^°Ca, -0.059 > > -0.069 fm and ^ospb 0.164 > 5„ > 0.074 fm. This encourages 
us to extend our global procedure to nuclei that have at least two data sets in the medium 
energy range. 

The sensitivity of Sn and Rn to the three different charge distributions, three parameter 
fermi (3PF), sum of gaussian (SOG) and fourier-bessel (FB) obtained from Ref. 6^ which 
are used in the fitting procedure is shown in Fig. Eland given in detail in Tables I, II and III. 
The five momentum transfer ranges, the proton form factor Gi{q) and the set FAOO Arndt 
NN amplitudes are used. For ^osp]^ g^^^^ ^^Ca there is little difference in the values of Sn and 
Rn for the different charge densities. We note that the RMS radii for ^''^Pb and ^^Ca for these 
three charge density are about the same jfi^l- Two of the three '^''Ca charge distributions 
also have almost identical charge RMS radii, 3PF 3.482(25)fm , SOG and 3.479(3)fm but 
the RMS radius for the FB, 3.450(10)fm, is considerably smaller (see Ref. jH^)- The result 
is that for ^°Ca the difference in the value of Sn is pronounced. As shown in Fig. El the FB 
charge distribution has range of Sn values that goes from positive to negative but the range 
of the Sn values 3PF and SOG are all negative. We attribute this to the smaller RMS radius 
for the FB charge distribution. The values of i?p, Rn and S'„ should not depend heavily 
on the momentum transfer ranges, the set of amplitudes, the form factors, and the charge 
distribution. In fact, the only case that does not overlap is the FB case for ^°Ca. The results 
for all three of the charge distributions are given in Tables I-III in the Appendix. 

In this global analysis we obtain many densities from the various cases used in checking 
the results of the fits. Figures 6 and 7 show the vector point proton density distribution, 

and r'^pP and the vector point neutron density distribution, and r^p" for "^^Ca, "^^Ca, 
and ^°^Pb. The SOG charge distribution and the Gi{q) form factor were used to obtain the 
vector point proton density. The neutron density also depends on the charge distribution 
chosen and the form factor used to obtain the vector point proton density as well as the other 
input to a given global fit, the NN amplitudes, and the EFT density used to fix the ratio 
of In Figs. 6 and 7 the Arndt NN amplitude set FAOO, the SOG charge distribution, 

the momentum transfer range was 0.0 fm~^ to 3.0 fm~^, proton from factor Gi(g), and all 
three EFT cases, MA4, FZ4 and VA3 are used. It is clear that the neutron densities overlap 
showing that the use of different EFT cases does not have any significant impact on the 
global fit. Tables of the densities and the parameters are available from the authors jfi^. 
This same input as in Figs 6-7 is also used in the Figs. 8-9 which shows the vector proton 
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FIG. 1: The ranges, shown as bars, of the skin thickness for ^'^Ca, ^^Ca and 208p]-, fj-Qj-j^ ^j^g 
global fits using the SOG charge distribution, set FAOO Arndt NN amplitudes, proton form factor 
Gi{q), and the three EFT models, MA4, FZ4 and VA3 are shown by the filled boxes, for the four 
momentum transfer ranges and the filled diamonds for the five momentum transfer ranges. Several 
theoretical values from Refs. 0, i g 0, H E3, 0, S, Q are also shown by using various 
symbols. 
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FIG. 2: The skin thickness with different proton form factors. Results for '^'^Ca, ^^Ca and ^''^Pb 
from the global fitting procedure shows the AVE eft models values for the SOG charge distri- 
bution, set FAOO Arndt NN amplitudes, the five momentum transfer ranges. We compare Gi{q) 
from Ref. [g^] shown as filled boxes and G2{q) from Ref. 61] shown as filled diamonds. The same 
theoretical values for Sn shown in Fig. ^ are also shown. 
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FIG. 3: The skin thickness for ^'^Ca, '^^Ca and ^'^'^Pb due to different sets of Arndt NN amphtudes. 
The global results shown are the AVEeft models for SOG charge distribution, the five momentum 
transfer ranges, and proton form factor Gi{q). Set FAOO is shown as filled boxes and set SM86 is 
shown as filled diamonds. The same theoretical values for Sn shown in Fig. ^ are also shown. 
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FIG. 7: The vector point proton density distribution shown by the sohd curve. The three neutron 
density distribution densities using the EFT densities, MA4, FZ4 and VA3 are shown by the dashed 
Une for case MA4, the dots for case FZ4, and the dotdashed for case VA3. All density are multiplied 
by the radius squared and obtained from the global fit as discussed in the manuscript. 
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point density and the scalar proton pf and r^pf, and in Figs. 10-11 for the vector point 
density and the scalar neutron p" and r^p". Figures 8-11 show that the phenomenology 
gives sensible scalar densities that are not unphysical. 

Figures 12 and 13 show the observables for one of the energies included in the global 
fit for ^°Ca and a prediction of an energy not included in the global fit. Figure 12 shows 
the fit for ^°Ca at 497.5 MeV and Fig. 13 shows the prediction for the ^°Ca 650 MeV spin 
observables; no cross section data is available. Figures 14 and 15 show the observables for 
one of the energies included in the global fit for ^^Ca and the results from a prediction of 
an energy not included in the global fit. Figure 14 shows the fit for ^^Ca at 497.5 MeV 
and the prediction of Q at that energy. Figure 15 shows the prediction for ^'^Ca at 1044 
MeV. Figures 16 and 17 show the observables for one of the energies included in the global 
fit for ^°^Pb and the results from a prediction of an energy not included in the global fit. 
Figure 16 shows the fit for ^osp]^ 497.5 MeV and Fig. 17 shows the prediction for the 650 
MeV observables for this target. The observables shown in Figs. 12-17 all used the Arndt 
NN amplitude set FAOO, the SOG charge distribution, the proton from factor Gi{q), the 
momentum transfer range from 0.0 fm"^ to 3.0 fm-\ and the EFT case MA4. 

While the figures of the observable are small we have magnified them and we find that 
the heights of the diffractive maxima and the angular positions of the minima and maxima 
are very well reproduced at each energy with no systematic energy dependent discrepancies. 
These are the most critical features of the data which determine rms radii. In fact, precision 
fits to Ay and Q are not as important for determining neutron radii but our fits are quite 
reasonable. The magnified figures are available from the authors j^. 

We also predict the total neutron cross section and the proton reaction cross section. In 
Fig. ^1 the predicted total neutron cross sections for ^°Ca and ^osp]^ jg compared with the 
experimental values from Finlay et al. [H^. The predicted proton reaction cross sections 
for the same two targets are given in Fig. with the experimental values j6^ and the 
predictions are reasonable. 

All of the results confirm that the global approach produces good fits to the data, that 
predictions are quite acceptable and that the neutron densities are energy independent. 
These densities should be useful input to a large number of nuclear reactions. They provide 
the empirical values of Rn and Sn that are needed for testing theoretical models. 



IV. SUMMARY AND CONCLUSIONS 

This paper reports a new method for extracting neutron densities from intermediate en- 
ergy elastic proton-nucleus scattering observables. Neutron densities are needed for atomic 
parity violation experiments, the analysis of antiproton atoms, the experiments parity vi- 
olation elastic electron scattering, and theoretical nonrelativistic and covariant mean-field 
models. It is interesting to note that our results for Sn are in agreement with work using 
antiprotonic atoms, see Fig. 4 in Ref. jl3]. 

The approach uses a global analysis, similar to the global fits using a Dirac phenomenol- 
ogy, but in this case it is based on the relativistic impulse approximation (RIA). The input 
to the procedure are the vector point proton density distribution, p*^, which is determined 
from the empirical charge density after unfolding the proton form factor and the Arndt NN 
amplitudes. The other densities, p", pf, p", are parametrized. The neutron densities for 
^"^Ca, ^^Ca and ^osp]^ obtained from the global fits have been tested for sensitivity to the 
input. Different proton form factors, different charge densities, different sets of the Arndt 
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FIG. 8: The vector point proton density distribution shown by the sohd curve. The three scalar 
proton density distribution densities using the EFT densities, MA4, FZ4 and VA3 are shown by 
the dashed hne for case MA4, the dots for case FZ4, and the dotdashed for case VA3. All density 
are obtained from the global fit as discussed in the manuscript. 
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FIG. 13: Prediction for ^^Ca at 650 MeV. 
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FIG. 14: Results of the global analysis for ^^Ca at 497.5 MeV. 
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NN amplitudes, different EFT models, different data sets included in the fits, and differences 
in the ranges of the momentum transfer used in the fits have been investigated. With the 
exception of the fourier-bessel charge density for ^°Ca all of these tests produce values for 
Rn and Sn that overlap. The prediction of data sets not used in the fits are well reproduced 
and the calculated observables are in good agreement with data. 

In conclusion, we have obtained values of Rn and Sn for ^"Ca, "^^Ca, and ^°^Pb. We have 
obtain good global fits for MA4, FZ4 and VA3 for every charge distribution used; SOG, 
3PF, and FB. Both five momentum transfer ranges and four momentum transfer ranges 
were used. As discussed earlier, for each case we calculate the mean and standard deviation 
for Rn and Sn for every global fit, this gives us an average over all EFT models for a 
given charge distribution. This is denoted by AVEeft in the tables. As mentioned several 
times this is not to be taken as a statistical analysis, but as a way to show the range of 
the values in a consistent way, they are shown as bars in the figures. We do this rather 
than remove outliers. Then for each case we show the combined values of the range for 
Rn and Sn for that case, this gives the average over all charge distributions as well as over 
all EFT models. These results are denoted as AVEcd and AVEeft in the tables. We 
then calculate the ranges of i?„ and Sn for all three cases combined for each target. From 
these final test cases, shown at the bottom of the tables in the appendix, using the five 
momentum transfer ranges we obtain the following values: for ^°Ca, 3.350 > Rn > 3.300 fm 
and -0.008 > Sn > -0.080 fm; for ^^Ca, 3.505 > Rn > 3.421 fm and 0.148 > Sn > 0.058 
fm; for ^ospb, 5.589 > Rn > 5.513 and 0.156 > S„ > 0.076 fm. In conclusion, the RMS 
neutron radius, Rn and the neutron skin thickness Sn obtained from the global fits that 
we consider to be the most reliable ranges of our results are given as follows: for ^°Ca, 
3.314 > Rn> 3.310 fm and -0.063 > ^„ > -0.067 fm; for ^^Ca, 3.459 > Rn > 3.413 fm 
and 0.102 > 5„ > 0.056 fm; and for ^o^Pb 5.550 > i?„ > 5.522 fm and 0.111 > Sn > 0.083 
fm. 

The authors have found that the global fits are more stable when using the Arndt NN 
amplitude set FAOO rather than the Arndt NN amplitude set SM86. This might be expected 
as the FAOO analysis is more recent. Both of the form factors used, Gi{q) and G2{q) when 
using FAOO were the most stable, see Tables I, II and III in the Appendix. The impact on 
the fits due to the charge distribution input, SOG, 3PF and FB, as shown in Fig. El clearly 
showed that the cases using the SOG and 3PF charge distributions were stable, i.e. the 
ranges overlapped, for all targets, however, for ^°Ca the results for the FB charge distribution 
gave very different results, see Table I. We also found that the fits using Arndt NN amplitude 
set FAOO, the Gi{q) form factor, and the SOG charge distribution produced the most stable 
results. In this case the values of i?„ and Sn using the five momentum transfer ranges are: for 
^°Ca, 3.318 > Rn> 3.308 fm and -0.059 > Sn > -0.069 fm; for ^^Ca, 3.498 > i?„ > 3.412 
fm and 0.141 > > 0.055 fm; and for ^o^Pb 5.602 > Rn > 5.512 fm and 0.164 > ^„ > 0.074 
fm. These values, as well as the values using the most conservative errors and the results that 
we consider most reliable, are in agreement with nonrelativistic Skyrme models jll li^. lisj]. 
relativistic Hartree-Bogoliubov model extended to include density dependent meson-nucleon 
couplings 0| and results from a recent analysis of antiprotonic atoms fl3\. Our results for 
^^Ca and ^osp^ ^re generally not in agreement with relativistic mean-field models, see Ref. 6 
and references there in. 

The global fit based on the RIA is a new tool for obtaining the neutron density. The 
values of Rn and Sn obtained are robust. This quality is verified in several ways as discussed 
in Section III. For example, we checked the sensitivity due to the data sets included in the 
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fit by doing global fits using only two data sets (497.5 MeV and 797.5 MeV) for ^°Ca and 
^°^Pb as well as four data sets. The values of Sn for the two data set case are: for '^''Ca 
-0.054 > ^„ > -0.066 fm; and for ^ospb 0.166 > Sn > 0.088 fm. These results agree 
well with Sn for the four data set case: for ^°Ca -0.059 > S„ > -0.069 fm; and for ^ospb 
0.164 > Sn > 0.074 fm. These results, shown in Fig. |1[ motivates us to use the global 
procedure for all nuclei that have at least two data sets in the medium energy range. 

This work provides energy independent values for i?„ and Sn, in contrast to the energy 
dependent values obtained by previous studies. In addition, the results presented in paper 
show that the expected rms neutron radius and skin thickness for ^°Ca is accurately repro- 
duced. The values of Rn and Sn obtained from the global fits that we consider to be the 
most reliable are given as follows: We plan to extend our work to additional nuclei and will 
continue to investigate different models and procedures. The goal is to continue to improve 
quality of neutron densities that result from our global fits. 
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TABLE I: 



Range (fm ^) 


Rn (fm) 


Sn (fm) 




«Ca Case 1 






Gl{q) 


FAOO 




SOG 


A^V^ pj' 


2.0-3.5 


3.314 > Rn> 3.310 


-0.063 > S„ > -0.067 


1.5-3.5 


3.318 > Rn> 3.308 


-0.059 > 5„ > -0.069 




3PF 


AVEeft 


2.0-3.5 


3.328 > R„> 3.304 


-0.051 > S„ > -0.075 


1.5-3.5 


3.328 > R„> 3.304 


-0.051 > S„> -0.075 




FB 




2.0-3.5 


3.373 > Rn> 3.349 


0.028 > S„ > 0.004 


1.5-3.5 


3.371 > Rn> 3.347 


0.025 > S„ > 0.001 




AVEcD 


-A V£/^jj'7^ 


2.0-3.5 


3.354 > R„> 3.306 


0.002 > 5„ > -0.076 


1.5-3.5 


3.352 > R„> 3.306 


0.000 > Sn> -0.076 


"'Ca Case 2 




G2{q) 


FAOO 




SOG 




2.0-3.5 


3.314 > iJ„ > 3.300 


-0.070 > S„ > -0.084 


1.5-3.5 


3.315 > Rn> 3.299 


-0.069 > Sn> -0.085 




3PF 


Ep p'j- 


2.0-3.5 


3.317 > Rn> 3.299 


-0.070 > S„ > -0.088 


1.5-3.5 


3.324 > i?„ > 3.296 


-0.062 > S„ > -0.090 




FB 




2.0-3.5 


3.341 > Ra> 3.333 


-0.012 > Sn> -0.020 


1.5-3.5 


3.340 > i?„ > 3.332 


-0.013 > > -0.021 






AVEepx 


2.0-3.5 


3.333 > R„> 3.301 


-0.027 > 5„ > -0.087 


1.5-3.5 


3.334 > ii„ > 3.302 


-0.026 > 5„ > -0.086 


40 Ca Case 3 






SM86 




SOG 




2.0-3.5 


3.344 > R„ > 3.308 


-0.033 > S„ > -0.069 


1.5-3.5 


3.341 > Rn> 3.305 


-0.036 > S„ > -0.072 




3PF 


A.V E^P'P 


2.0-3.5 


3.332 > Rn> 3.302 


-0.047 > 5„ > -0.077 


1.5-3.5 


3.330 > R„ > 3.292 


-0.049 > Sn> -0.087 




FB 




2.0-3.5 


3.374 > R„> 3.340 


0.028 > 5„ > -0.006 


1.5-3.5 


3.371 > Rn> 3.335 


0.025 > Sn> -0.011 








2.0-3.5 


3.356 > ii„ > 3.308 


0.002 > Sn> -0.070 


1.5-3.5 


3.354 > Rn> 3.304 


-0.001 > 5„ > -0.075 


^Ca for all three cases 


2.0-3.5 


3.350 > Rn> 3.304 


-0.006 > S„ > -0.080 


1.5-3.5 


3.350 > iin > 3.300 


-0.008 > Sn> -0.080 
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TABLE II: 


Range (fm ^) 


Rn (fm) 


Sn (fm) 




48 Ca Case 1 






Gl{q) 


FAOO 




SOG 


pj' 


2.0-3.5 


3.459 > Rn> 3.413 


0.102 > S„ > 0.056 


1.5-3.5 


3.498 > Rn> 3.412 


0.141 > Sn> 0.055 




3PF 




2.0-3.5 


3.456 > -R„ > 3.412 


0.087 > S„ > 0.043 


1.5-3.5 


3.501 > Rn> 3.409 


0.132 > S„> 0.040 




FB 




2.0-3.5 


3.460 > -Rn > 3.420 


0.113 > S„ > 0.073 


1.5-3.5 


3.499 > Rn> 3.417 


0.152 > S„ > 0.070 




AVEcD 




2.0-3.5 


3.459 > Rn> 3.415 


0.103 > Sn> 0.055 


1.5-3.5 


3.499 > Rn> 3.413 


0.142 > Sn> 0.054 


^^(. 'ii. (.'aso 2 




G2{g) 


FAOO 




SOG 


AVEeft 


2.0-3.5 


3.453 > -R„ > 3.407 


0.089 > S„ > 0.043 


1.5-3.5 


3.498 > Rn> 3.404 


0.134 > Sn> 0.040 




3PF 


AVEp P'j- 


2.0-3.5 


3.448 > R„ > 3.404 


0.072 > S„ > 0.028 


1.5-3.5 


3.500 > Rn> 3.398 


0.124 > S„ > 0.022 




FB 


AVEepx 


2.0-3.5 


3.458 > il„ > 3.414 


0.103 > 5„ > 0.059 


1.5-3.5 


3.511 > R„> 3.407 


0.157 > S„ > 0.053 




AVEcD 


AVEepx 


2.0-3.5 


3.454 > Rn> 3.408 


0.092 > 5„ > 0.040 


1.5-3.5 


3.503 > Rn> 3.403 


0.139 > 5„ > 0.037 


*^Ca Case 3 




Gi(</) 


SM86 




SOG 


A\^ E E FT 


2.0-3.5 


3.506 > Rn> 3.476 


0.150 > S„ > 0.120 


1.5-3.5 


3.505 > > 3.463 


0.148 > S„ > 0.106 




3PF 


A.V E^P'P 


2.0-3.5 


3.487 > Rn> 3.447 


0.119 > 5„ > 0.079 


1.5-3.5 


3.484 > Rn> 3.446 


0.116 > Sn > 0.078 




FB 


AVEeft 


2.0-3.5 


3.516 > -R„ > 3.486 


0.169 > S„> 0.139 


1.5-3.5 


3.515 > Rn> 3.473 


0.167 > S„ > 0.125 






AVEeft 


2.0-3.5 


3.509 > Rn> 3.465 


0.158 > 5„ > 0.100 


1.5-3.5 


3.504 > il„ > 3.458 


0.152 > Sn> 0.094 


**Ca for all three cases 


2.0-3.5 


3.485 > -R„ > 3.417 


0.129 > S„ > 0.053 


1.5-3.5 


3.505 > Rn> 3.421 


0.148 > Sn> 0.058 
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TABLE III: 



Range (fm ^) 


Rn (fm) 


Sn (fm) 




208P6 Case 1 






Gl{q) 


FAOO 




SOG 




2.0-3.5 


5.550 > R„> 5.522 


0.111 > S„> 0.083 


1.5-3.5 


5.602 > R„ > 5.512 


0.164 > 5„ > 0.074 




3PF 


AVEeft 


2.0-3.5 


5.580 > -R„ > 5.548 


0.155 > > 0.123 


1.5-3.5 


5.586 > Rn> 5.552 


0.160 > Sn> 0.126 




FB 




2.0-3.5 


5.546 > -Rn > 5.518 


0.108 > Sn> 0.080 


1.5-3.5 


5.600 > Rn> 5.508 


0.162 > 5„ > 0.070 




AVEcD 




2.0-3.5 


5.565 > Rn> 5.523 


0.135 > 5„ > 0.085 


1.5-3.5 


5.599 > Rn> 5.521 


0.167 > Sn> 0.085 


-"^rij Case 2 




G2(?) 


FAOO 




SOG 


AV Ee ft 


2.0-3.5 


5.545 > -R„ > 5.517 


0.102 > 5n > 0.074 


1.5-3.5 


5.599 > Rn> 5.507 


0.156 > 5„ > 0.064 




3PF 


AV Ep PJ' 


2.0-3.5 


5.576 > -R„ > 5.538 


0.146 > Sn > 0.108 


1.5-3.5 


5.582 > Rn> 5.542 


0.152 > Sn> 0.112 




FB 




2.0-3.5 


5.537 > P„ > 5.511 


0.095 > Sn> 0.069 


1.5-3.5 


5.593 > R„> 5.499 


0.151 > Sn > 0.057 




AVEcD 


AVEeft 


2.0-3.5 


5.558 > Rn> 5.516 


0.124 > Sn > 0.074 


1.5-3.5 


5.594 > Rn> 5.514 


0.157 > Sn > 0.075 


208P6 Case 3 




Gi(</) 


SM86 




SOG 


AVEe ft 


2.0-3.5 


5.538 > R„ > 5.498 


0.100 > Sn > 0.060 


1.5-3.5 


5.554 > Rn> 5.502 


0.115 > Sn > 0.063 




3PF 




2.0-3.5 


5.584 > Rn> 5.532 


0.158 > Sn > 0.106 


1.5-3.5 


5.590 > Rn> 5.538 


0.165 > Sn > 0.113 




FB 




2.0-3.5 


5.532 > Rn> 5.502 


0.094 > Sn > 0.064 


1.5-3.5 


5.551 > Rn> 5.503 


0.113 > Sn > 0.065 






AV Eeft 


2.0-3.5 


5.559 > Rn> 5.503 


0.129 > Sn > 0.065 


1.5-3.5 


5.571 > il„ > 5.509 


0.141 > Sn > 0.071 


^o^Pb for all three cases 


2.0-3.5 


5.561 > -R„ > 5.513 


0.130 > Sn > 0.074 


1.5-3.5 


5.589 > Rn> 5.513 


0.156 > Sn > 0.076 
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